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In a recent letter [Phys. Rev. Lett. 82, 1012 (1999)] it was found that the Tb'^"'" magnetic moments 
in the Tb2Ti2 07 pyrochlore lattice of corner-sharing tetrahedra remain in a collective paramagnetic 
state down to 70mK. In this paper we present results from d.c. magnetic susceptibility, specific heat 
data, inelastic neutron scattering measurements, and crystal field calculations that strongly suggest 
that (1) the Tb^+ ions in Tb2Ti2 07 possess a moment of approximatively 5/iB, and (2) the ground 
state tensor is extremely anisotropic below a temperature of O(10'')K, with Ising-like Tb^^ mag- 
netic moments confined to point along a local cubic (111) diagonal (e.g. towards the middle of the 
tetrahedron). Such a very large easy-axis Ising like anisotropy along a (111) direction dramatically 
reduces the frustration otherwise present in a Heisenberg pyrochlore antiferromagnet. The results 
presented herein underpin the conceptual difficulty in understanding the microscopic mechanism(s) 
responsible for Tb2Ti2 07 failing to develop long-range order at a temperature of the order of the 
paramagnetic Curie-Weiss temperature 9cw ~ — lO^K. We suggest that dipolar interactions and 
extra perturbative exchange coupling(s) beyond nearest-neighbors may be responsible for the lack 
of ordering of Tb2Ti2 07. 
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I. INTRODUCTION 

In magnetic systems, competition between magnetic in- 
teractions, combined with certain local lattice symme- 
tries involving triangles, give rise to the notion of geo- 
metric frustration . Geometrically frustrated antifer- 
romagnets are currently attractingmuch interest within 
the condensed matter community . The main reason 
for this interest is that geometric frustration can cause 
sufficiently large zero-temperature quantum spin fluctu- 
ations as to drive a system into novel types of intrinsi- 
cally quantum mechanical magnetic ground states with 
no classical equivalent |||-|^ . 

Among three-dimensional systems, the pyrochlore lat- 
tice of corner-sharing tetrahedra (see Fig. 1) with an- 
tiferromagnetic nearest-neighbor exchange interaction is 
particularly interesting. For this system, theory 
and Monte Carlo simulations show that for classical 
Heisenberg magnetic moments interacting with a nearest- 
neighbor antiferromagnetic coupling, there is no transi- 
tion to long-range magnetic order at finite temperature. 
This is unlike the two-dimensional kagome lattice anti- 
ferromagnet |^,|io|,0 where a thermally-driven order-by- 
disorder of spin nematic order occurs. Villain coined the 
name "collective paramagnetic" to describe the classical 
state of the pyrochlore lattice at low temperatures H . 




FIG. 1. Pyrochlore lattice of corner-sharing tetrahedra. 
The magnetic moments occupy the corners of the tetrahedra. 

Because of their low propensity to order even for clas- 
sical spins, antiferromagnetic materials based on a py- 
rochlore lattice appear to be excellent systems in which 
to seek exotic quantum mechanical ground states. For 
example, numerical calculations suggest that the S = 
1/2 pyrochlore Heisenberg antiferromagnet may be fully 
quantum disordered, giving rise to a state that is com- 
monly referred to as "spin liquid" Both terms 
"collective-paramagnet" and "spin liquid" are meant to 
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emphasize that despite such a system remaining in a 
paramagnetic phase down to absolute zero temperature, 
the properties of such a state involve very strong and 
nontrivial short-range spin correlations, analogous to the 
nontrivial position-position correlations present in an or- 
dinary atomic or molecular fluid. 

A number of experimental studies on insulating py- 
rochlore materials have been reported in the past ten 
years. Interestingly, it has been found that such 
systems do not typically form such a spin liquid 
state that remains paramagnetic down to zero tem- 
perature. Most often, these systems either display 
long-range antiferromagnetic order; such as FeFs [12|, 
GdsTiaOr IHjll, ZnFc204 and ZnCr204 Sjig, 

or enter a spin-glass-like state below some nonzero spin- 
freezing temperature as exhibited by Y2M02O7 pS^EOj, 
Tb2Mo207 Bill, Y2Mn207jii, as weh as the disor- 
dered CsNiCrFg pyrochlore |2^. 

Recently, several studies of the pyrochlore rare-earth 
titanates, R2Ti207, have been published I2II. In these 



compounds, the trivalent rare earth ions, K , occupy 
the 16d sites of the FdSm space and form a pyrochlore 
lattice (Fig. 1). The behaviors displayed in this fam- 
ily of pyrochlores are much varied indeed. Gd2Ti2 07 
develops true long-range order at a critical temperature 
of about IK |ll3| , |l4| . Tm2Ti207 possesses a trivial non- 
magnetic (i.e. spin singlet) ground state separated by 
an energy gap of about 120 K to the next crystal field 
level ||25|1 . Ho2Ti207 is well described by an Ising dou- 
blet |2q| . In that system, it was originally argued that 
the nearest-neighbor exchange interaction is weakly fer- 
romagnetic , and that the strong Ising- like single ion 
anisotropy along (111) directions frustrates the develop- 
ment of long range ferromagnetic order |p6|^-|33|| . This ma- 
terial also exhibits low temperature spin dynamics remi- 
niscent of Pauling's "ice model" 13^ , an equivalence pro- 
posed by Harris and co-workers jl26| , ^ -|3^ . Recently, it 
has been found that Dy^TizOv is also a very good 
example of "spin ice" 1 26 , and that the applica- 

tion of a magnetic field can restore much of the ground- 
state entropy and drive magnetic phase transitions. Most 
recently, den Hertog and Gingras have argued that the 
spin ice physics in both IIo2Ti207 and Dy2Ti2 07 is not 
driven by nearest-neighbor ferromagnetic exchange, but 
is rather due to the long-range nature of magnetic 
dipole-dipole interactions p^ , pO[ |. 

In contrast to the long-range ordered or spin glass 
states mentioned above, strong evidence for collective 
paramagnetism, or spin liquid behavior, was recently ob- 
served in the insulating pyrochlore Tb2Ti207 [^6 37|. It 
was found using neutron scattering and muon spin re- 
laxation methods that this material remains paramag- 
netic down to (at least) 70 mK despite the fact that the 
paramagnetic Curie- Weiss temperature, ^cw, is -19 K, 
and that short-range antiferromagnetic correlations be- 
gin to develop at ~ 50 K. At first sight, one could argue 
that it is "pleasing" to have found at last the spin liq- 
uid state anticipated by theory for a highly frustrated 



pyrochlore antiferromagnet |^-^. However, the situation 
for Tb2Ti207 is not as simple as it might naively appear. 

In Tb2Ti207 the Th^+ ions have a partially filled '^Fg 
shell, and one most first understand their crystal field 
level scheme and, in particular, the nature of the single- 
ion magnetic ground state before constructing a correct 
effective spin-spin Hamiltonian for Tb2Ti207. Indeed, 
we show below that crystal field anisotropy renders the 
description of Tb2Ti207 in terms of an isotropic Heisen- 
berg antiferromagnetic model completely inappropriate. 

If one neglects the axial oxygen distortion around the 
Tb'^"'" sites, and assumes that the local environment of the 
Tb'^"'" is perfectly cubic, one would expect, based on point 
charge calculations, that the ground state of both Tb^"*" 
and Tm'^^ should either be a singlet or a nonmagnetic 
doublet p^ . For example, as mentioned above, experi- 
mental evidence for a nonmagnetic singlet ground state 
has been found in Tm2Ti207 (25|. Based on this naive 
picture, one can see that the experimental evidence of a 
moment for Tb'^+ in Tb2Ti207 is therefore a nontrivial 
issue that needs to be understood. 

A simple possibility is that corrections beyond the 
point-charge approximation and/or the known axial oxy- 
gen distortions around each of the 16d sites cause the 
Tb'^+ cations to acquire a permanent magnetic moment. 
Another and more interesting possibility, is that the mo- 
ment on the Tb'^^ site in Tb2Ti207 is induced by a col- 
lective bootstrapping of the magnetic (exchange and/or 
dipolar) interactions as occurs in the tetragonal LiTbF4 
material 1^^. In Tb2Ti207, a priori, it is theoretically 
possible that there could be no moment on the Tb'^^ 
site for a concentration x of Tb^+ less than some criti- 
cal concentration Xc in (Tb2;Yi_a;)2Ti207, as occurs in 
LiTbj;Yi_a:F4 This is an important issue. In- 

deed, one could imagine that for the highly-frustrated 
pyrochlore lattice, the collective development of a per- 
manent ground-state moment would not give rise to ho- 
mogeneous moments on the Tb"^^ sites, but to a kind of 
"modulated moment structure" . This idea is conceptu- 
ally similar to what is found in the frustrated tetragonal 
TbRu2Ge2 material but where for Tb2Ti207 there 
might be instead a quantum-disordered state "interven- 
ing" between a trivial singlet ground state and a long- 
range ordered one, with the quantum-disordered state 
extending all the way to x — 1. 

In the case where a permanent moment does exist on 
Tb'^^ even in absence of interaction (i.e. the limit x ^ 
in (Tba;Yi_a;)2Ti207), the important issue is to deter- 
mine the wavefunction decomposition of the ground state 
in terms of | J, AIj) states and the symmetry, Heisenberg 
or otherwise, of the resulting effective spin variable. The 
goal of such a programme is to construct a low-energy ef- 
fective spin Hamiltonian in order to tackle theoretically 
why Tb2Ti2 07 does not order at nonzero temperature. 
Consequently, it is very important to know in more de- 
tails what the magnetic nature of the Tb'^^ single-ion 
ground state in Tb2Ti207 is. 

The main purpose of this paper is to examine the mag- 
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netic nature of the Tb^+ ion in the Tb2Ti207 pyrochlore 
in order to assess whether or not there is indeed a perma- 
nent moment at the Tb site as the temperature goes to 
zero, and determine the nature of this moment (e.g. level 
of effective spin anisotropy) . We present in Section II ex- 
perimental evidence, based on results from d.c. suscepti- 
bility, heat capacity and powder inelastic neutron studies 
that show there is a permanent moment at the Tb site, 
but that its approximate S/xb value is less than the value 
of 9.4/.tB estimated from the d.c. susceptibility measure- 
ments above 200K or the 9.72/xb Fq free ion value. 
To complement the experimental work, results from ab- 
initio crystal field calculations that take into account co- 
valent and electrostatic effects are presented in Section 
III and Appendix B. We discuss in Section IV the pos- 
sibility that dipole-dipole interactions and extra pertur- 
bative exchange couplings beyond nearest-neighbor may 
be responsible for tlie lack of ordering in Tb2Ti207. 



II. EXPERIMENTAL METHOD & RESULTS 
A. Sample Preparation 

Samples of Tb2Ti207 and (Tbo.o2Yo.98)2Ti207 were 
prepared in the form of polycrystalline pellets by high 
temperature solid state reaction. Starting materials, 
Tb2 03, Y2O3 and Ti02 were taken in stoichiometric 
proportions, mixed thoroughly, pressed into pellets and 
heated in an alumina crucible at 1400°C for 12 hours 
in air. Tb2 03 was prepared by hydrogen reduction of 
Tb407. The powder x-ray diffraction patterns of the 
samples, obtained with a Guinier-Hagg camera, indi- 
cate that they are single phase with cubic unit cell 
constants, ao, of 10.49lA for Tb2Ti207 and 10.1041 
for (Tbo.o2Yo.98)2Ti207. The value for the concen- 
trated sample is in excellent agreement with previous 
reports @42|. 

Some of this high quality polycrystalline material was 
then used as starting material for a successful single crys- 
tal growth using an optical floating zone image furnace. 
Details of the crystal growth are described elsewhere [Esl . 



B. d.c. Magnetic Susceptibility Measurements 

As a first step towards determining the magnetic na- 
ture of the electronic ground state of the Tb^+ cations in 
Tb2Ti207, we have investigated the d.c. magnetic sus- 
ceptibility of Tb2Ti207 and (Tbo.o2Yo.98)2Ti207. The 
d.c. magnetic susceptibility was measured using a 
SQUID magnetometer(Quantum Design, San Diego) in 
the temperature range 2-300 K. The inverse susceptibil- 
ity, x~^5 of Tb2Ti207 measured at an applied field of 
0.01 Tesla is shown in Fig. 2. A fit of the data to 
the Curie- Weiss (CW) law above 200 K gives an effec- 
tive paramagnetic moment of 9.4/iB/Tb'^+ and an effec- 



tive (high-temperature) CW temperature, ^cw = -19 K, 
which indicates the dominance of antiferromagnetic in- 
teractions (the Curie- Weiss fits are done above 200 K). 
Deviation from the Curie- Weiss law sets in at a rather 
high temperature, T ~70 K. This is consistent with clas- 
tic neutron scattering results reported previously, which 
showed evidence for short range magnetic correlations at 
temperatures up to 50 K |3q| . 

Recognizing that for Tl?+, an ^Fg even electron and 
non-S-state ion, there may be "crystal field" as well as 
exchange contributions to the experimentally determined 
6'cw, a- magnetically dilute sample, (Tbo.o2Yo.98)2Ti207, 
was also studied. The data for this sample are shown 
in Fig. 3 along with a Curie- Weiss fit giving again a 
value close to the free ion value for the effective param- 
agnetic moment, 9.6/^B/Tb'^+, and 6'cw ~ —6 K. This 
finite value is in contrast to the essentially zero ^cw 
value obtained for a similarly diluted sample of Gd2Ti207 
which contains the isotropic "spin only" ^S7/2 Gd^+ ion 
]l3[ , and therefore indicates that a significant crystal field 
contribution to 6'cw exists in the Tb^+— based material. 
Thus, to a first approximation one can estimate that the 
portion of ^cw for Tb2Ti207 which can be attributed to 
magnetic interactions is 

0Cw{Tb2Ti2O7 } - ecw{(Tbo.02Yo.98)2Ti207 } ^ -13 K, 

a value similar to the 6'cw ~ -10 K found for 
Gd2Ti207 This approach can be made more rig- 

orous by noting that in a high-temperature series ex- 
pansion, one finds that the magnetic susceptibility is 
X = Ci(l/T C2/T2) where C2 = 6'cw can be "de- 
composed" as a simple sum of terms that are ascribed to 
exchange interactions, dipolar interactions and crystal- 
field terms (see Appendix A). Note also that, down to a 
temperature of T = 2 K, neither Fig. 2 nor Fig. 3 show 
any sign of a singlet ground state which would manifest 
itself as an approach to a constant susceptibility with 
decreasing temperature, as found in the single ground- 
state of the Tm2Ti207 |25|. A more detailed analysis 
concerning this issue is presented in Section III. 

While it is tempting to use the magnetic interaction 
contribution (^cw — ^cw) ~ — 13 K, (with 6'cw ~ — 19 
K and 6'^vv ~ — 6 K), to extract the approximate value of 
the nearest neighbor exchange, an estimate of the near- 
est neighbor dipole-dipole interaction for Tb'^+ ions in- 
dicates an energy scale of about 1 K. Given the long 
range nature of dipolar forces, it becomes clear that the 
classical nearest neighbor exchange constant cannot be 
obtained from a measurement of the Curie- Weiss tem- 
perature until the effects of long range dipolar interac- 
tions on ^cw have been understood. We find via a high 
temperature series expansion analysis of the long range 
dipolar contributions, 6'^^, to the Curie- Weiss tempera- 
ture (see Appendix A), that the estimated upper bound 
on is ferromagnetic and ~ 1.2 K (for needle-shaped 
powder crystallites), while the lower bound is antiferro- 
magnetic and ^ —2.4 K (for slab-shaped powder crystal- 
lites). Consequently, we find that antiferromagnetic ex- 
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change interactions are predominantly responsible for the 
{9cw - Ocyf) = -13 K value determined above T > 200 
K, with a resulting 6'cw ^ [-14.2, -10.6] K. 
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FIG. 2. Inverse molar susceptibility, 1/x, of Tb2Ti2 07 vs 
temperature. 
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FIG. 3. Inverse molar susceptibility, of the diluted 

(Tbo.o2Yo.98)2Ti207 material vs temperature. 



C. Neutron Scattering Experiments 

Inelastic neutron scattering allows us to determine 
with reasonable precision the values of the electronic en- 
ergy levels of Tb2Ti207. The inelastic neutron scatter- 
ing measurements were carried out on a 50g sample of 
polycrystaUine Tb2Ti207, loaded in a sealed Al cell with 



a helium exchange gas present. The cell was mounted 
in a closed cycle helium refrigerator with a base tem- 
perature of 12 K. The sample was the same as used in 
reference p6| . Measurements were performed on the C5 
triple axis spectrometer at the Chalk River Laboratories 
in constant scattered neutron energy mode. Two spec- 
trometer configurations were employed, appropriate for 
relatively high and low energy resolution, respectively. 
Both configurations employed pyrolitic graphite (PG) 
as both monochromator and analyser. The low resolu- 
tion measurements, appropriate for relatively high energy 
transfers, were performed using -^=3.52 THz (lTHz=48 
K), open-60-80-open collimation, and a PG filter in the 
scattered beam. The high resolution configuration used 
^=1.2 THz, open-40-60-open collimation, and a cooled 
Be filter in the scattered beam. These results clearly in- 
dicate excitations at i? w 0.35THz, 2.5THz and a broad 
neutron group centered at 3.5THz (corrsponding to 16.8 
K, 120 K, and 168 K, respectively). 

The low energy-resolution inelastic neutron scattering 
measurements at 12 K with -^=3.52 THz revealed the 
presence of two Q-independent modes with frequencies 

~ 2.5 and 3.5 THz. Representative neutron groups, as 
well as the dispersion relation for these two excitations 
are shown in the top panel of Fig. 4. These excitations 
are identified as being magnetic in origin due to their 
temperature and Q dependence. Their flat dispersion in- 
dicates that they are crystal electric field levels for Tb'^+ 
in the environment appropriate for Tb2Ti207. 

The high energy-resolution inelastic neutron measure- 
ments at 12 K with ^-—1.2 THz shows the presence 
of a low lying magnetic excitation near v ~ 0.35 THz. 
This mode is also dispersionless above a temperature of 
~ 25 K, but partially softens in energy at the wavevec- 
tor which characterizes the very short range spin cor- 
relations, which develop below 25 K. The development 
of this interesting dispersion has been described previ- 
ously Js^ . The dispersion of the low lying excitation is 
shown in the lower panel of Fig. 4. One can see that this 
partial softening of the excitation branch occurs only for 
the lowest lying mode. As we show in the next section, 
this is manifested in the heat capacity measurements as 
broad features that result of a broadening of the single- 
ion energy levels via these magnetic correlation effects. 

These measurements place constraints on any calcula- 
tions for the energy eigenstates of Tb'^"'" as they set both 
the energy spacing of the levels, and require that mag- 
netic dipole matrix elements must connect the ground 
state with these levels in order that they be visible to 
the inelastic neutron scattering experiment. In particu- 
lar, this indicates nonzero (0| J+|l), (0| J-|l), or (0| J^|l) 
matrix elements connecting the ground state, |0), and 
the first excited state, |1), at an energy ~ 0.35 THz « 
17 K above the ground state. In other words, there must 
be large \J,Mj) components in |0) and |1) where some 
of the Mj involved for the ground state and the excited 
state differ by 0, ±1. 
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FIG. 4. Inelastic neutron scattering data obtained at a 
temperature T = 12 K. The inset of the top panel of Fig. 4 
shows a transfer energy scan at Q = 2A~^ . Data showing rep- 
resentative neutron groups as well as the dispersion relation 
for these excitations at ~ 2.5 and 3.5 THz is shown in the top 
panel of Fig. 4. The dispersion of the low lying excitation is 
shown in the lower panel. 
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FIG. 5. Specific heat, Cp, of Tb2Ti207 as function of tem- 
perature for the indicated magnetic fields of 0, 2 and 5 T. The 
solid line corresponds to the lattice specific heat of Tb2Ti2 07, 
Ci, estimated from the measurements on the non- magnetic 
Y2Ti207 that is isotructural to Tb2Ti207 [52]. 



D. Specific Heat Measurements 

Low-temperature specific-heat measurements on 
Tb2Ti2 07 were performed using a thermal- relaxation mi- 
crocalorimeter. The single crystal sample was mounted 
on a sapphire holder which was isolated from the bath by 
four copper-gold alloy wires. The relative precision and 
absolute accuracy of the calorimeter were confirmed by 
measuring copper and gold standards. In principle, spe- 
cific heat-measurements on a dilute (Tbo.o2Yo.98)2Ti207 
sample would be useful. Unfortunately, this is not tech- 
nically easily feasible as the magnetic contribution to the 
total specific heat would be too small to be determined 
accurately. 



The total specific heat, Cp, of Tb2Ti207 was measured 
from 0.4 K to 30 K at appHed fields of 0, 2 and 5 T (see 
Fig. 5). The zero field data exhibits two broad peaks cen- 
tered at about 1.5 K and 6 K. The data are in agreement 
with those of Ref. |^^, above a temperature of 0.4 K. 
Hyperfine contributions to the specific heat become im- 
portant below 0.4 K for Tb— based compounds, as found 
for example in the Tb2(GaSn)07 pyrochlore [^^, and 
this is presumably the reason for the sharp increase of 
Cp{T)/T found in Fig. 4 of Ref. (27) below 0.5 K. The 
solid line in Fig. 5 corresponds to the estimated lattice 
heat capacity, C/, for Tb2Ti207, determined by scaHng 
the heat capacity for Y2Ti207, which is insulating, non- 
magnetic, and is isotructural to Tb2Ti207 [ p2[ . 

The magnetic specific heat, C„i, obtained by subtract- 
ing Ci from Cp, is shown in Fig. 6 for the three applied 
fields. With the application of a 2 T field the magnitude 
of the lower temperature peak diminishes and moves to a 
slightly higher temperature. In contrast, the position of 
the second peak does not change but increases in magni- 
tude as the lower temperature feature begins to overlap 
with it. At 5 T the low temperature feature disappears 
completely and the remaining peak is shifted to higher 
temperatures. 
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FIG. 6. Magnetic specific heat, Cm, of Tb2Ti207 as func- 
tion of temperature for the indicated magnetic fields of 0, 2 
and 5 T. The solid lines are guides to the eye. 

The crystal field calculations described in Section III 
and Appendix B (see Table 2) indicate a level scheme 
consisting of a ground state doublet with another doublet 
as the first excited state. Attempts to fit the Cm data 
to Schottky anomalies using a unique doublet— doublet 
level scheme for the Tb'^^ ions and varying the ground 
state— excited state energy splitting failed, since impor- 
tant magnetic short-range correlations are present in this 
sytem, as discussed above in neutron scattering [ p6[ . 
Consequently, we interpret the anomaly at ~ 6 K as 
a remnant of the excited doublet that is "broadened" 
by exchange correlation fields, while the 1.5K anomaly 
is presumably due to these same correlation effects, but 
now acting on the single-ion ground state doublet. In 
other words, the build-up of short-range correlations in 
the low-temperature sector of the effective Hamiltonian 
for Tb2Ti207 results in a specific heat anomaly at 1.5 K. 
This low-temperature anomaly at 1.5K that results from 
correlations is akin to the broad specific heat bump at 
« 2 K found in GdaTiaOr (iS). However, GdaTiaOy is 
an spin-only ion, and there are no crystal field lev- 

els at high energy, nor are there correlation remnants of 
crystal field levels above the ground state such as those 
that cause the specific heat anomaly at 6 K in Tb2Ti207. 
Another distinction between Tb2Ti207 and Gd2Ti207 is 
that the latter shows a very sharp specific heat anomaly 
at 0.9 K |l^, and therefore presumably a transition to 
long range order at that temperature as suggested by 
recent neutron scattering experiments [^. Down to 
0.4 K, no such sharp specific heat anomaly is found in 
Tb2Ti207. From the fit of the d.c. susceptibility and the 
crystal field calculations presented in Section III and Ap- 
pendix B, good evidence is obtained that the magnetic 
moment in the ground state and the first excited state 



is ~ 5 /zb and ~ 6 /xb, respectively. Given a doublet- 
doublet energy gap of about 17 K, we can estimate the 
strength of the magnetic field where the separation be- 
tween the doublets is equal to the magnetic field energy, 
and find a magnetic field of about 5 T. For an applied 
field of that strength, the ground state and excited states 
merge and are strongly mixed. This explains the disap- 
pearance of the low-temperature specific heat anomaly 
in Fig. 6 for a fuekd H — 5 T. 

It is usual to estimate the magnetic entropy, Sm{T), 
in a system by integrating Cm{T)/T between the lowest 
temperature reached and the temperature, T, of interest. 
Although it is straightforward to integrate Cm{T)/T, the 
interpretation of the results for Tb2Ti207 is difficult. 
The main reasons for this are: 

• The hyperfine interaction is large for Tb'^"'" and the 
nuclear specific heat contribution becomes signifi- 
cant with respect to the magnetic contribution near 
0.4K |0|j5|]. 

• There is a doublet crystal field excitation at an en- 
ergy of approximately 17 K. Hence, one can hardly 
integrate Cm{T)/T above 5 K without "already" 
embedding in the resulting entropy a contribution 
from excitations to the first excited doublet. 

• From our neutron results and crystal field calcu- 
lations presented in Appendix B, it is known that 
there are other levels at an energy ~ 10^ K. Hence 
one cannot integrate Cm(r)/T to obtain Sm{T) 
up to a high enough temperature without having 
to consider the specific heat contribution from the 
states at - 0{W) K. 

• By integrating Cm{T)/T up to ~ 30 K one enters 
a regime where the lattice contribution to the spe- 
cific heat, Ci{T), becomes sizeable (see Fig. 5). In 
that case the subtraction of Ci{T) from the total 
Cp{T) using rescaled results for the isostructural 
non- magnetic Y2Ti2 07 leads to inherent uncertain- 
ties which increase dramatically above 10 K (see 
Fig. 5). 

With these provisions in mind, we have determined 
Sm{T) (see Fig. 7). We find that the recovered entropy, 
Sm{T), at 15 K is already larger than that expected for a 
singlet— doublet energy level scheme, S'„i(15K) > Rln(3). 
Since 15 K is much less than the excited states at ~ 100 
K, there should be little contribution to ^^(ISK) com- 
ing from states at T > 100 K. Consequently, the results 
presented in this figure support further the picture that 
the two lowest energy levels in Tb2Ti2 07 consists of two 
doublets. However, it is interesting to note that at a tem- 
perature of 30 K, the recovered entropy is not yet equal to 
Rln(4), the total entropy for two doublets. Either 30 K is 
not yet at high enough temperature to have recovered the 
full doublet— doublet entropy, or there exists macroscopic 
entropy in the ground state as occurs in Dy2Ti207 
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For the four reasons mentioned above, it is difficult to 
make the discussion about the recovered entropy above 
30 K in Tb2Ti207 much more quantitative. 

In summary the magnetic specific heat data are consis- 
tent with the inelastic neutron scattering results in that 
the two lowest-lying energy levels for Tb^+ in Tb2Ti207 
consist of two doublet energy levels separated by an ex- 
citation energy of ~ 15 — 20 K. As we will show in the 
next Section and Appendix B, our crystal field calcula- 
tions strongly suggest that these two lowest lying levels 
are magnetic (Ising) doublets. 
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FIG. 7. Magnetic specific heat, Cm(T), of Tb2Ti207 in 
zero applied field (circles), and entropy, Sm{T) (squares) are 
shown. The entropy values for two, three and four states are 
indicated by the horizontal dashed lines. R is the molar gas 
constant R= Noks =8.3145 J mol~^ K~^. The solid lines are 
guide to the eye. 



III. SINGLE ION PROPERTIES 

As discussed in the introduction, it is very important 

to determine whether the existence of a moment at the 
Tb'^+ site in Tb2Ti207 is intrinsic, or driven by mag- 
netic interactions. Consequently, we have investigated in 
detail the problem of single- ion properties of Tb^+ , both 
theoretically and via d.c. susceptibility measurements, of 
Tb2Ti207 and of the dilute (Tbo.o2Yo.98)2Ti207 mate- 
rial where the Tb— Tb interactions should be negligible. 
The main conclusion from this investigation is that the 
Tb^"*" cation does indeed carry an intrinsic moment in- 
herent to its environment in Tb2Ti207, and that moment 
is not due to a bootstrapping effect from interactions. 

The results from inelastic neutron scattering and spe- 
cific heat measurements presented in the previous sec- 
tion already provide some evidence for a doublet ground 



state and an excited doublet state at an energy of ~ 
17K. The purpose of this section is to investigate further 
from a theoretical point of view the question of the ex- 
istence of a magnetic doublet ground state for Tb'^+ in 
Tb2Ti207. In short, both a simple point charge calcula- 
tion and a more sophisticated ab-initio method confirm 
that a doublet— doublet scheme is the most likely low en- 
ergy level structure for Tb^"*" . We have constructed a van 
Vlcck equation based on such a doublet— doublet scheme 
in order to parametrize the d.c. suceptibility of the dilute 
(Tbo.o2Yo.98)2Ti207 sample, and to extract the value of 
the magnetic moment for both the ground state and ex- 
cited state doublets. We find that experimental results 
for the d.c. susceptibility of (Tbo.o2Yo.98)2Ti207 are in 
reasonably good agreement with the ab-initio calcula- 
tions. 



A. Crystal Field Effects 

Experimental evidence of a doublet-doublet structure 
for the low temperature crystal field levels of Tb^"*" in 

Tb2Ti207 can be understood by considering the crystal 
field environment surrounding the ion. Pyrochlore ox- 
ides A2B2O7 are described in space group Fd3m with 
A^+, the trivalent rare earth in 16d, B'*+, the tetravalent 
transition metal ion, in 16c, 01 in 48e and 02 in 8b. The 
A or Tb^"*" site in this c;asc is coordinated to six 01 ions 
at about 2.5A in the form of a puckered ring and to two 
02 ions at a distance of 2.2A in the form of a linear 02- 
Tb-02 chain oriented normal to the mean plane of the 
01 ring. The 02-Tb-02 units are parallel to the (111) 
directions within the cubic unit cell. Overall, the local 
geometry at the Tb'^+ site can be described as a severe 
trigonal compression along the body diagonal of a simple 
cube. 

Based on symmetry considerations, the cubic plus axial 
distortion surrounding the Tb'^^ ion may be expressed by 
a crystal field Hamiltonian of the general form 



fT-cf _ r2^2 , r4^4 , h4,^4 
+ -Dq^O + -03^3 + -06^6 



(3.1) 



where the B^'s are yet to be determined crystal field 
parameters, and the Cg 's are tensorial operators defined 
as Cg = [47r/(2fc + 1)]^/^Y^, where is a normalized 
spherical harmonic. In general, the B^'s represent an 
effective one-body potential which lifts the degeneracy of 
the angular momentum states in question. In practice, 
they may be determined experimentally by spectroscopic 
or thermodynamic probes, or theoretically within various 
levels of approximation. 

In the simplest approximation of the crystal field in- 
teractions, one often uses the so-called point charge (PC) 
approximation where the crystal field is simply assumed 
to be caused by the Coulomb field of point charges situ- 
ated at neighboring sites. In such a picture, the crystal 
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field eigenstates and eigenvalues have been determined 
for a number of rare earths by Lea et al. for systems with 
cubic symmetry (i.e. without trigonal distortion) 
For Tb^+, the lowest three energy levels in a cubic en- 
vironment are a Tt, singlet, a non-magnetic T2 doublet 
(2) . 

and a Fg triplet, with their precise ordering in terms 
of lowest energy states dependent upon variation of the 
point charge crystal field parameters Pq] . In general, the 
addition of the trigonal distortion splits the Fg ' triplet 
into a singlet and doublet, while the F3 and F2 states are 
preserved. Hence, we expect that the crystal field ground 
state of the Th'^^ ion in Tb2Ti207 will be a competition 
between two doublets and two singlets. 

Some notion regarding the difficulty of determining the 
precise ordering of these states and the size of their as- 
sociated magnetic moments may be obtained by diago- 
nalizing the crystal field Hamiltonian iJ^^ using Stevens' 
operator equivalents of the within a fixed J man- 
ifold [Q, and by using the point charge approximation 
for the Coulomb effects of the surrounding oxygen ions. 
The resulting crystal field point charge Hamiltonian, Hp^^, 
with quantization axis along the appropriate (111) direc- 
tion can be expressed as, 



H;i = ajBlir' - 1)0° 



ijBt 



(27r^ 



1: 



28 

(188-h324r^ 
512 



^6 



35V2 



ol 



77 



^6 



(3.2) 



where r — R1/R2 and i?i, i?2 are the Tb-0 distances for 
oxygen ions situated on the puckered ring and on the 
distortion axis respectively. The O"^ represent crystal 
field operators as discussed by Hutchin gs |5q| , while aj, 
and jj are the Steven's coefficients [ f33y55[ ] . Trivalent 
Tb^"*" is an '^Fq ion and thus the fixed J manifold is J = 6 
{L = 3, S = 3, J = L + S = 6) foi the operator equiv- 
alent point charge calculation. The precise relationship 
between the point charge parameter set {B^} and the 
more general {B^} is discussed by Kassman | |54| . 

Although a point charge estimation of the crystal field 
parameters is in most cases unreliable in predicting the 
actual crystal field level spacing of rare earth ions, we 
find that varying the point charge parameter set {B^} 
indicates that the ground state can be confirmed to be 
a competition between two singlets close in energy and 
two magnetic doublets, which are well separated from 

the other crystal field states. These levels are indeed 

(2) 

the remnants of the F2,F3 and Fg states of the cubic 
environment eluded to earlier. 

In general, we find that in a large region of the crystal 
field parameter space, the two doublets form the low- 
est energy levels, although their precise ordering may 
change. Although the weight of each angular momentum 
component of a crystal field eigenstate varies with the 
values of the crystal field parameters, some other general 



features do emerge. In particular, one singlet contains 
only I ± 6) and | ± 3) states while the other is a combina- 
tion of I ± 6) , I ± 3) and 1 0) . On the other hand both sets of 
doublets are magnetic (e.g they have a nonzero quantum 
expectation value of J^), and also contain components 
of exclusively different values (the operator does 
not connect the ground state to the excited state). One 
doublet has large | ± 4) and | ± 1) components, while 
the other has large | ± 5) and | ± 2) components. Hence, 
the two doublets have AIj components that differ by ±1, 
and a neutron spin-flip induced transition from one to 
the other is allowed, consistent with what is found in the 
inelastic neutron scattering results presented above. 

We have confirmed these conclusions based on our 
point charge analysis by performing a more sophisticated 
first principles calculation that take into account both 
electrostatic and covalency effects as well as the intra 
atomic and configurational interactions. This approach, 
described in Appendix B, does not restrict the decompo- 
sition of the electronic energy levels into a fixed \J,Mj) 
manifold as is usually done using the Steven's operator 
equivalents, as discussed above. In the results presented 
in Appendix B, we find that the two lowest energy dou- 
blets have a leading Mj = ±4 and Mj ± 5 components, 
respectively (close to 90% of the weight). Table 2 in Ap- 
pendix B lists three very similar energy level structures 
given for slightly different constraints on the crystal field 
parameters. The theoretically determined energy levels 
(mean values: 0, 13, 60 and 83 cm~^, that is 0, 19, 86, 
and 119 K, respectively, be compared with the experi- 
mental levels determined by inelastic neutron diffraction 
(0, 17, 115 and 168 K). The corresponding wavefunctions 
allow for a good estimate of the magnetic susceptibility 
of the diluted and concentrated compounds as shown by 
Table 3 and Figs. 2, 3, and 8. 

From these results, a picture of the low temperature 
single-ion properties of Tb^+ can be deduced. Consid- 
ering the structure of the eigenstates of the two lowest 
doublets, a calculation of their 5— tensors indicates ex- 
tremely strong Ising like anisotropy along the appropri- 
ate (111) axes (the axis formed by joining the two cen- 
ters of the tetrahedra that the ion belongs to) for each 
Tb'^"'" ion at low temperatures. In summary, based on 
both point charge and ab inition calculations, we have 
strong evidence of a doublet— doublet scheme for Tb^^^, 
and a substantial single-ion anisotropy in Tb2Ti207 p^ ] 
making the the Tb^+ moment effectively Ising— like for 
T < 0(10"'^) K. The consequences of this result will be 
discussed in Section IV. 



B. Susceptibility 

The zero field susceptibility measurements of the pow- 
der Tb-diluted compound (Tbo.o2Yo.98)2Ti207 shown in 
Fig. 3 enables us to gain a more concise understanding of 
the nature of the crystal field levels of Tb'^"'' in Tb2Ti207, 
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and the single ion properties of Tb^+. Indeed, the dilute 
concentration of Tb'^'^ions in the system removes the ef- 
fect of magnetic interactions and in principle leaves only 
crystal field contributions to the magnetic susceptibility. 
Having obtained strong evidence for a doublet— doublet 
scheme for the Tb^^ ions we now proceed further and 
analyse the d.c. magnetic susceptibility by construct- 
ing a phenomenological expression for the susceptibility 
based on a van Vleck equation for such a doublet— dou- 
blet energy level structure (7|]. This allows us to ex- 
tract the size of the magnetic moments of both doublets 
from experiment, as well as check for consistency with 
the conclusions based on our crystal field calculations of 
Appendix B. 

Due to the powder nature of our (Tbo.o2Yo.98)2Ti207 
sample, the random orientation of the grains can lead to 
both parallel and transverse contributions to the suscep- 
tibility at first and second order. For the doublet-doublet 
structure at low temperature eluded to earlier, the van 
Vleck equation for the susceptibility has the general form 

ffVliVTb ( a/T + b + e-^/^{c/T d) \ 



1 



where g is the Lande factor, equal to 3/2 for Tb'^"'', /is 
is the Bohr magneton and TVxb = 0.04iVoj where Nq is 
Avogadro's constant. The adjustable parameters a, b, c, d 
are defined through second order perturbation theory to 
be 



\mo) 



= E l(«il^"l'^i> 



where nQ,mo label states within the ground state dou- 
blet. Til, mi label states in the excited doublet while the 
index i defines any state from the ith crystal field level. 
The A's represent crystal field energy level differences 
(in Kelvin) and a = x,y or z. The fitting parameters a 
and c are due to first order terms in perturbation theory 
while b and d are from second order terms and give rise 
to temperature independent van Vleck paramagnetism 
contributions to x- 

We have performed a least squares fit to the suscepti- 
bility data of (Tbo.o2Yo.98)2Ti207 up to a temperature 
~ 30 K, which is approximately where thermal contribu- 
tions from a crystal field level at ^ 100 K become non- 
negligible. Because of the narrow energy spacing between 
the ground state and excited state doublet, all four ad- 
justable parameters a, b, c, d are important in fitting the 
susceptibility data. Based on the crystal field calcula- 
tions of the previous section, the specific heat analysis of 
the doublet-doublet gap in the concentrated Tb2Ti2 07 
sample, and the evidence of an anisotropy gap of 17 K 
observed in inelastic neutron measurements on the same 
sample |36 , we have carried out the fit of the low tem- 



perature behavior of the susceptibility using a doublet- 
doublet gap A ranging from [12 ^ 24] K. We find that 



the goodness of fit is quite flat in this range for A. How- 
ever the magnitude of the adjustable parameters do not 
deviate strongly and can be determined to a reasonable 
degree of accuracy over this interval. Using values of 
the gap outside of this interval yields a noticeably poorer 
goodness of fit. In Fig. 8, we show the fit to susceptibil- 
ity data for the dilute (Tbo.o2Yo.98)2Ti207 sample using 
an anisotropy gap of 17 K as well as the results for the 
ab-initio crystal field calculations. 

We can interpret the fitted results for the susceptibility 
data by making use of our crystal field results and Eq. 
( |3.3[ ). Due to the strong Ising like nature of the g— tensors 
of the two theoretically calculated doublets, transverse 
terms such as {ni\J^\mi) between two states within a 
doublet are negligible compared to {ni\J^\ni). Thus, 
we expect that a ~ |('t-o|'^^|^o)P and c ~ |(ni|J^|ni)p 
and consequently, both the a and c terms represent per- 
manent moment contributions to the susceptibility. For 
the same reason, the t;- tensor characterizing the ground 
state is extremely anisotropic with essentially only a 
component along the local (111) direction with very little 
g± component. As a result, and this is the most impor- 
tant point of the paper: 

At a temperature T < 0(10^) K, Tb^+ ions can be 
considered to a very good approximation as (classical) 
Ising magnetic moments confined to point parallel or 
antiparallel to their local (111) directions. 

Over the range mentioned above for the doublet-doublet 
gap A, we have found that the magnitude of the moment 
in the doublet ground state to be 5.10 ± 0.3^b. Overall, 
we find that our best fit for the susceptibility data gives 
a = 11.6 ±.1,6= 1.53 ± .04 K-\c = 15.7 ± 4.0 and 
d = .71 ± .05 K^^. The value of c gives a magnitude for 
the moment of the first excited doublet of 5.9 ± .S/i^. 

Both fitted moments are compatible with the eigen- 
state structures of the doublets determined from the crys- 
tal field calculations presented in Appendix B. Addition- 
ally, the values of the paramagnetic terms b and d in the 
susceptibility are also consistent with our crystal field re- 
sults. In particular, from the J'^ components of our cal- 
culated low energy doublet and singlet eigenstates, there 
will be predominant contributions to b and d coming from 
transverse angular momentum matrix elements connect- 
ing the two doublets, as well as additional contributions 
involving transverse matrix elements between the dou- 
blets and the higher energy singlet states at > 100 K. 
The doublet-doublet coupling will give equal (in mag- 
nitude) contributions to b and d while coupling to the 
singlets will give further positive contributions to b while 
reducing the value of d. This can be simply understood 
in terms of the signs of the denominators for each of these 
virtual excitation processes in our definitions of b and d. 

In summary, we are able to successfully fit the 
(Tbo.o2Yo.98)2Ti207 susceptibility measurements at low 
temperature using a doublet-doublet picture consistent 
with our crystal field calculations. We find reasonable 
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agreement in terms of calculated moments and param- 
agnetic contributions between that expected from theory 
and our fitted values. In particular, we find a perma- 
nent moment in the ground state of approximately 5.1/iB- 
This moment is intrinsic to the Th'^^ ion and is not 
driven by exchange and/or dipolar interactions as occurs 
in LiTbF4 This value is also compatible to what is 
estimated from the limiting low-temperature muon spin 
relaxation rate 1 /Ti found in Ref. |M] , assuming a dipole 
coupling between a positive muon bounded to an oxy- 
gen at ~ 2.5A away from a ~ 5/iB Tb'^"'" moment (see 
Appendix C). 
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FIG. 8. Inverse susceptibility of (Tbo.o2Yo.98)2Ti207. Cir- 
cles in dica te measured values while the solid line is a fit using 
Eq. (3.3) using a doublet-doublet gap of A = 18 K. The 



sohd triangles are the values obtained from the theoretical 
calculations of Appendix B. 



IV. DISCUSSION 

Combining together our results from d.c. susceptibility 
measurements, specific heat data, inelastic neutron scat- 
tering data, and crystal field calculations, the following 
picture emerges: 

• The Tb'^+ ion in Tb2Ti207 carries a permanent 
magnetic moment of approximately 5/iB ■ That mo- 
ment is intrinsic to the Tb'^"'" ion and is not driven 
by magnetic correlation from exchange and/or 
dipolar interactions. 

• The ground state is well described as an Ising dou- 
blet with extremely anisotropic (7— tensor. In other 
words, the moments in the single-ion ground state 
are predominantly confined to point along the local 
(111) directions. 

• The first excited state is at an energy of approxi- 
mately 18 ±1 K above the ground state, which is 



also characterized by a very anisotropic Ising-like 
g— tensor. 

• Consequently, in the absence of interactions, the 
Tb^+ ions should be very well modeled below T < 
O(IO^) K as effective classical Ising spins confined 
to point along the local (111) directions. 

From d.c. susceptibility measurements at high temper- 
ature, we know that the magnetic interactions are pre- 



iCxchangc 
CW 



-13 K. 



dominantly antiferromagnetic, with 6\ 
Although the interactions are not small compared to the 
first excitation energy gap of 18 K, let us momentarily 
ignore the exchange coupling between the ground state 
and the excited doublet, and consider only the Ising-like 
ground state doublet. We then have a classical model 
with effective Ising spins pointing along their local (111) 
directions. In other words, spins that can only point ei- 
ther inward or outward on a tetrahedron. 

Because of the open structure of the pyrochlore lat- 
tice, we expect that the nearest-neighbor exchange in- 
teractions predominate. If we make the further approx- 
imation that only nearest-neighbors contribute to the 
antiferromagnetic interactions, we arrive at a scenario 
where Tb2Ti207 is effectively described by classical Ising 
spins pointing along the (111) directions and coupled 
via nearest-neighbor antiferromagnetic exchange. As dis- 
cussed in Refs. pl|j7^ , [7^ , such a model is not frustrated. 
Indeed, if we pick one spin, and chose either an "in" or 
"out" orientation for it, then the three other spins on 
the same tetrahedron must also take the same "in" or 
"out" configuration to minimize the exchange interac- 
tions. Since the pyrochlore lattice can be described as an 
FCC lattice with a tetrahedron as the basis cell, either 
an "all in" or "all out" configuration repeats identically 
on all basis units; one refers to such a magnetically long- 
range ordered structure as a Q = ground state |f7|,pO[. 

The observation that Tb2Ti207 remains paramagnetic 
down to at least 70 mK, while sustaining important short 
range magnetic correlations, is therefore very puzzling. 
We believe that the most plausible explanation is that the 
interactions in Tb2Ti2 07 involve further than nearest- 
neighbor interactions, J„„. The presence of these interac- 
tions, as well as the long-range dipolar interactions, pos- 
sibly "reintroduce" large frustration in the spin Hamil- 
tonian describing Tb2Ti207, and conspire to destroy the 
"would-be" long-range Neel Q = ground state, and give 
rise to a collective paramagnet— spin liquid, ground state. 



V. CONCLUSION 

In conclusion, we have presented results from d.c. sus- 
ceptibility, specific heat, inelastic neutron scattering, and 
crystal field calculations for the pyrochlore lattice antifer- 
romagnet Tb2Ti207. We have obtained strong evidence 
that the Tb'^^ magnetic moment on the 16d site at T < 2 
K is intrinsic and is not induced by magnetic (exchange 
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and/or dipolar) interactions or correlation effects such as 
found in LiTbF4 |3§], and is of the order of 5/xb ■ AH ev- 
idence points towards a very strong Ising like anisotropy 
for the doublet ground state which forces the resulting 
classical Tb'^^ Ising moments to point either parallel or 
antiparallel to their local (111) direction. For antiferro- 
magnetic exchange interactions, such strong anisotropy 
largely removes all local ground state spin degeneracy, 
and should naively force the system to possess an Ising- 
like long range order ground state with all spins in or out 
on a tetrahedron basis cell. The reason for the failure 
of TbaTisOy to order at a temperature ~ [10° - 10^] K 
set by the exchange part of the Curie- Weiss temperature 
remains unresolved at this time. The results presented 
here point towards the need to consider the role that ex- 
change interactions beyond nearest-neighbor and dipolar 
interactions play in Tb2Ti207. 
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APPENDIX A: 
Dipole Contributions to ^cw 

The contribution from dipolar interactions to ^cw niay 
be determined from a high temperature series expansion 
of the long range dipole-dipole Hamiltonian to first order 
in (3. That is to say, given the additive property of crystal 
field effects, exchange and dipolar interactions to 9cw, 

i.e. 9cw = ^'cw 
from an expansion of 



and consequently that 



gex 



^CW' may determine 



2 2 

kT 



where z defines some global direction and we assume T 3> 
dew isotropic spins. The angular brackets reflect 
an expectation value with respect to the dipole-dipole 
Hamiltonian iJ'^'P defined by 



^dip 



1 



2, ,2 \ " 
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To first order in /3 this yields 
g^filNJiJ + l 



X 



3kT 



^ gV|j(J+l) ^^ 
3NkT 



ndip 

where 



gV|j(J+l) 
3iVfc 



A 



(A.l) 
(A.2) 



We see from the above analysis that the evaluation of 
involves the summation of a conditionally conver- 
gent series (a lattice sum), and thus special care must be 
taken in its evaluation. In general, it can be treated in 
a controlled manner by the use of a rapid convergence 
factor via the Ewald method @|5) . Within the Ewald 
approach, the sum is split into two rapidly converging 
sums, one over the real space lattice and one in recip- 
rocal space. Additionally, a surface (shape dependent) 
term also arises (which is interpreted as a demagnetiza- 
tion factor ||,|6|). 

Indeed, if we approximate the powder grains of the 
(Tbo.o2Yo.98)2Ti207 as spherical, then the sum Eq. (A.2) 
is identically zero. This can be understood by imagining 
Eq. (A.2) as a sum over parallel (along the z axis) dipoles 
moments of magnitude unity. For a spherical object, it 
can be shown that such a sum must be identically zero 
for a system with cubic symmetry 1 47 4^ . Accordingly 
there would be no dipolar contribution to the measured 
value of dew |49|. On the other hand, we do not expect 
that the geometry of the powder grains is in fact spher- 
ical, (additionally there is also the possibility of effects 
from inter-granular interactions). Therefore to gain an 
estimate on the upper bound of the dipolar contribution, 
we carry out the lattice sum Eq. (A.2) for an infinitely 
long cylinder (needle shape) along the z direction where 
surface effects are zero. This allows us to gain an approx- 
imate upper bound on the dipolar contribution as there 
are effectively no demagnetization effects. This calcula- 
tion can be carried out rather simply by noting that for 
a spherical sample, Eq. (A.2), can be written using the 
Ewald method as [Esl : 



Aspherc 
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L bulk 



A bulk I A surface 
^'-sphere 
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M ^ 
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The Rij(n) = r^/L + n, i,j label the dipole (Tb''+) 
sites within the cubic unit cell, L is the length of the 
conventional cubic cell for Tb2Ti207 and M is the 
number of cells in the sample. The function H{y) — 
(2ly^pK) dx while n = {k,l,m) such that k,l,m 

are integers and a is suitably chosen to ensure rapid con- 
vergence of the summations over n. Based on the pre- 
viously mentioned result that Agpher should vanish for 
a spherical sample, this implies that the first four terms 
in the above expression, A''""^, must sum to the oppo- 
site value of the last term (A^p^^"'^'^, the surface term), 
and indeed we have verified that this is the case numer- 
ically. For an infinitely cylindrical sample, the surface 
term, A^^j'j^^^J: = 0, and thus one finds 



A A 

-^^-cylindcr 



_ \bulk _i_ ^surface 



cylinder 



^bulk 



M ^ 47r 

M 

-1005— r . 



(A.6) 



On the other hand, for a slab geometry, we have 

A A bulk I A surface 

Aslab = A -I- A„ 



slab 



(A.7) 



where 
thus arriving at 

Aslab 



A surface ^ \ ^ ^ 
slab =;^Z^4^ 



M ^ Stt 
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Combining Eqs.(A.l), (A.6) and (A. 8), and using a cubic 
cell length of L 10.104A, we arrive at upper and 
lower bounds for dipolar contributions to the Curie- Weiss 
temperature, namely —2.4 < 6*^'^ < +1.2 K. 



APPENDIX B: 
Crystal Field Calculations 

Our aim in this Appendix is to pursue in more details 
and using a more sophisticated (ab-initio) approach the 
calculations of the magnetic susceptibilites of Tb2Ti207 
and (Tbo.o2Yo.98)2Ti207. To do so, we need to deter- 
mine the electronic structure of the ground ^Fg level of 
the Tb^"'' (4/* configuration) in these two compounds, 
deduce the wavefunctions and from there infer the value 
of the magnetic moment of the ground state level. The 
determination of the 4/^ electronic configuration is ob- 
tained by diagonalizing the following Hamiltonian for 
a generic /" configuration without making the fixed J 
manifold approximation used in Section III. A: 

HiD ^ J2 F'iff)f' + C(/)^o + c.L{L + 1) + 



kq 



(B.l) 



Let us explain first what are the different terms in _ff (/"). 

• The i^'^'s [k = 2, 4, 6) are the electrostatic integrals 
(Slater's parameters) which splits the 4/" configu- 
rations into terms where S is the total spin 
and L is the total orbital angular momentum. The 
/*'''s are the associated two-electron operators 

• C(/) is the spin-orbit interaction integral which 
splits the terms into ^'^+^Lj levels. Aso is the as- 
sociated one-electron spin-orbit operator |^6| . 

• a, (3 and 7 are parameters associated with effective 
two-body correction terms for inter-configuration 
interaction |5^. C(G2) and C(R7) are the 
Casimir's operators for groups G2 and R7. When 
2 < n < 12, there can be several terms ^^'^^L with 
the same S and L values in the f" configuration. 
For instance there are three terms in 4f® while 
there exists only one term ^F. The states may dif- 
fer by the way they are built from the parent con- 
figuration f"^^. An additional classification of the 
states is therefore necessary. It is done according to 
the irreducible representations of the the groups G2 
and R7 and bestows additional quantum numbers 
to the states. 



the are the coefficients of the one-electron 
crystal field interaction which acts between 
P'^+^iJMj) sublevels. They can be theoretically 
predicted or extracted from fits of the energy lev- 
els (spectral lines) from experiments. In the point 
charge electrostatic model, their expression is : 



i?^^ = (4V2fc +1)1/2 (,'=)i]^.(Q^./i?- 



Where 



is a 4/ electron radial integral, Qj is 



the point charge of ligand j, and Rj, Oj, and are 
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the polar coordinates of ligand j. The derivation of 
the for the covalent interactions is much more 
involved The = (47r/2k+l)i/2[[Y«]], are 

the tensorial onc-clcctron crystal field operators. 

The evaluation of the matrix elements of i) the elec- 
trostatic interaction, ii) the spin-orbit interaction, Hi) 
the free-ion configuration interaction, and iv) the crys- 
tal field interaction, between the states of the basis set 
chosen for the /" configurations is necessary in order to 
determine the eigenvalues and eigenvectors of the lat- 
ter. The matrix elements are calculated by the means 
of tensorial algebra |Q. Besides, if coupled p'^+^LJA/j) 
states are chosen as the basis set, the one- or two-electron 
operators which are involved in the Hamiltonian cannot 
act directly on them. The calculation requires interme- 
diary mathematical quantities known as reduced matrix 
elements which are tabulated for standard /" configura- 
tions ( 1^). Once evaluated, the complete matrix ele- 
ments are multiplied by the associated parameters before 
diagonalization. The parameters are then determined by 
trial and error by successive diagonalizations and com- 
parison of the eigenvalues with experimental energy lev- 
els. In the present case our main concern is the structure 
of the Tb^+ ground level, so that the only specific mate- 
rial dependent parameters which have to be determined 
before diagonalizing are the B^ crystal-field in- 

teraction parameters. As already mentioned , the most 
convenient way to deduce the electronic structure of a 
rare earth ion in a solid compound, and hence determine 
the crystal-field parameters, is usually via analysis of its 
electronic spectrum by fitting the B^ to match the fre- 
quency of the observed transitions. The strongest lines 
are due to electronic transitions partly allowed by the 
mixing of the ground configuration with opposite par- 
ity configurations ( Judd-Ofelt mechanism) [^^jIiI • How- 
ever, in the case of pyrochlores, the site symmetry at 
the rare earth site is centrosymmetrical. The odd par- 
ity crystal field parameters vanish, the mixing of oppo- 
site parity configurations is impossible, hence no electric 
dipole transitions are detectable in the spectrum. In- 
deed, only the weak ^Fi^^Dq and ^Fq^^Di magnetic 
dipole transitions were observed previously in the absorp- 
tion spectra of the pyrochlore compounds Eu2Ti207 and 
Eu2Sn207 Hence a complete set of "phenomeno- 

logical" crystal field parameters (CFP) cannot be deter- 
mined from optical absorption or emission spectra. How- 
ever, recent inelastic neutron scattering experiments have 
been able to give information on the lowest electronic lev- 
els of Ho2Ti207 From these neutron results, some 
CFP's can be deduced. CFP's can also be calculated ab- 
initio from the compound structure and the atomic data 
of the constituents. Therefore, in what follows, two ap- 
proaches are used to determine the CFP's of Tb2Ti2 07 
and (Tbo.o2Yo.98)2Ti207: 

I A) A full ab-initio calculation of the CFP's uti- 
lizing the structural data of the compounds. 



IB) The fit of the CFP's from the Ho2Ti207 ^ 
and Eu2Ti207 data and a transposition 
to the Tb2Ti207 and (Tbo.o2Yo.98)2Ti207 
compounds. 

The next two steps for the calculations of the magnetic 
susceptibility are: 

II) The calculation of the 4/^ electronic configu- 
ration utilizing plausible free ion parameters, 
and the fitted or calculated ab-initio values of 
the CFP's. 

Ill) The calculation of the magnetic susceptibility 
of Tb2Ti207 and (Tbo.o2Yo.98)2Ti207 utiliz- 
ing the wavefunctions derived from the pre- 
vious step II. 

These steps are detailed in the following : 

lA) Ab-initio calculation of the CFPs 

An ab-initio determination of the CFP's is obtained by 
adding an electrostatic and a covalent contribution along 
lines similar to the ones developed in Refs. |^,^ for 
oxygen ligands. The crystal structure, the ionic charges 
and the ionization energies of the ligands are used. In 
Rcf. [|64| , experimental and predicted values of the pa- 
rameters calculated by the "covalo-electrostatic" model 
were compared for ten compounds. The mean deviation 
between experimental and calculated values: 



AB'^/B'"' 



— k<.q<.k —k<q<k 
[l + l/(5'^)2-2cOs(i?'=)/5'=]l/2 



1/2 



where and R}' are the scale and reliability factors 
listed in Table 7 of Reference Q. AB'^/B'' is found 
to be equal to 52, 30 and 23% for k ^ 2, k ~ 4 and 
k = 6, respectively. Such is the uncertainty which can be 
expected from a "blind eyed" prediction of the CFP's of 
Tb2Ti207. 

As mentioned earlier, the space group of rare earth ti- 
tanates with the pyrochlore structure is Fd3m. The eight 
oxygen first neighbours form a distorted cubic polyhe- 
dron. Two oxygens occupy ideal positions on opposite 
summits of the cubic threefold axis. The three [[ sides ]] 
of the cube originating from each of these two summits 
are equally elongated. The a cubic lattice parameter is 
equal to 10.15 and 10.09 A, and x the positional param- 
eter for the six displaced oxygens to 0.3 and 0.2968 for 
the dense and dilute compound respectively |42) . 

The site symmetry at the rare earth site is reduced 
from Oh (cubic) to D^d- The remaining threefold order 
symmetry axis imposes for the crystal field parameters 
the condition g = 0, modulo 3, so that the non-zero even 
k CFP's are B^, B^, B^, B^, B| and B|. The predicted 
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CFP values are reported in Table 1. The distances be- 
tween the Tb'^"'" ion and the oxygens on the threefold 
axis is short (2.20 and 2.18 A for the dense and the di- 
lute compound, respectively) while the distances to the 
six peripheral oxygens are much larger (2.52 and 2.49 A, 
respectively). This explains why the "axial" Bq param- 
eters are much larger than the "azimuthal" B^'s . 

IB) Experimental determination of the CFPs. 

The CFPs can be determined by fitting the energy lev- 
els of the /" configuration either using spectroscopic or 
inelastic neutron scattering data from other pyrochlore 
compounds. Here, we can use the Bq from spectroscopic 
data on Eu2Ti207 ||6|] and the other B^ from neutron 
scattering data on Ho2Ti207 ||27|] . 

• Determination of Bq from spectroscopic data on 
Eu2Ti207 

The ''Fi^^Do and '^Fq— »^Di magnetic dipole tran- 
sitions were observed previously in the electronic 
absorption spectrum of Eu2Ti2 07 In the lat- 
ter, the ^Fi and ^Di splittings amount to 291 K 
and 51.9 K, respectively. A fit in 4/^(Eu^+) yields 
Bq = 684 K. The transposition to Tb'^^ is made as- 
suming the crystal field parameter is scaled by the 
ratio of radial integrals: B^{Th^+) = B^{Eu^+) 
x(r2)(Tb3+)/(r2)(Eu3+) = Sg(Eu3+) x 0.91 = 
622 K. This "experimental" value, listed in Table 
1, and referred to as (a), has the same sign, but it 
is about half of the predicted value. As pointed out 
hereabove, the uncertainty on the calculated k — 2 
parameters is large. 

• CFP determination in Ho2Ti207. 

Siddharthan et al. recently reported results 
from ineleastic neutron scattering at low tempera- 
ture in Ho2Ti207. They determined the six lowest 
E irreducible representations of the ^Ig ground state 
level. Utilizing their experimental values, we fitted 
the CFP's of Ho^+ while maintaining the ratio be- 
tween CFP's with the same k value close to the 
theoretical ratio. The program ATOME was used 
for the refinement |^^. In this program, the basis 
set is composed of Slater determinants which makes 
unnecessary the use of tables of reduced matrix el- 
ements []65| . The evaluation of the matrix elements 
is straightforward, but the configuration cannot be 
truncated. Indeed, each eigenvector being a linear 
combination of a large number of Slater determi- 
nants, none of the latter can be omitted. As a con- 
sequence, all the 1001 states of 4f^° configuration 
are included in the diagonalization matrix. The ba- 
sis is large (1001) but still tractable. The final mean 
deviation between experimental and calculated lev- 
els was equal to 7.8 K. The fitted CFP's are : B^^ 
3173 K, B|= -1459 K, B§= 1343 K, Sf= 1292 



Kand i?|= 609 K. As pointed out hereabove, the 
CFP's were then scaled according to the ratio be- 
tween the radial integrals of Ho'^"'" and Tb^"*" to give 
the experimental CFP's for Tb2Ti207. Namely, 
B^(Tb3+) = B^-(Ho3+) x(r'=^)(Tb3+)/(r'=^)(Ho3+). 
These "experimental values" are are listed in Table 
1, and referred to as (b). The experimental k = 4 
and fc = 6 CFP's are 1.7 and 1.9 times larger than 
the theoretical values which is somewhat unusual. 

Ill) Calculation of 4/^ electronic configuration. 

The calculation of the 4/^ electronic configuration 
is done by the means of program /" utilizing 
the Hamiltonian H{f") in B.l acting on coupled states 
^^~^^Lj. Contrary to program ATOME previously men- 
tioned, program P can work on a truncated basis, which 
is necessary to resolve the 4/*configuration of Tb'^+with 
a large number of states (3003 states). The interac- 
tion matrix is built on a 387x387 basis set comprising 
the following ^■^+^L terms of the Tb'^+( 4/®) configura- 
tion: 7f, 5d(1,2,3), 5f(1,2), 5G(1,2,3), 3p(l,2,3,4,5,6) 
and ^8(1,2,3,4). The conjugate configuration of Tb^"*" is 
that of Eu3+ with 4/ + 2-8 = 6 electrons. The if (n=6) 
configuration of Eu^+ contains exactly the same number 
of basis states (e.g. [4/-|-2]!/[[4/-|-2-n]!n!] = 3003 states) 
as Tb'^^, the same terms, and the same levels. The in- 
teractions involving an odd number of electrons have a 
reverse sign in 4f^ and 4f®. For instance the ^■^+^Lj lev- 
els appear in a reverse order, and so do the crystal field 
sub-levels. In addition, the terms determined by the elec- 
trostatic interaction (two electron interaction) appear up 
in the same order for Eu^+ and Tb'^^. 

For Eu^+ the above quoted basis had proved large 
enough to allow a simulation of the levels up to ^D2 
(30219 K) without drastic truncation effects The 
F'^'s were assigned the Gd^+ values given in Ref. [|68| , 
that is 147289 K, [[ 102479 K ]], and 55868 K for fe = 
2, 4 and 6, respectively, a, (3, 7 were ascribed the Nd"^"*" 
values fitted in Ref. [|9|, that is 30.98 K, -1005.03 K, 
and 2510.48 K, respectively. The spin-orbit coupling con- 
stant (■(/) was set equal to 2446.30 K which is a stan- 
dard value for Tb'^+ Bq was assigned the trans- 
posed value quoted hereabove, and the other CFP's the 
values listed in Table 1 obtained after rescaling the B^ 
extracted from fits to the levels of Ho'^+. The diagonal- 
ization of the interaction matrix gives the energy levels 
and the corresponding leading eigenvectors listed in Ta- 
ble 2. The lowest levels are two doublets. Both states are 
rather Ising-like, with nearly exclusive (~ 0.95) Mj = ±4 
and Mj — ±5 components. This is caused by the very 
"axial" crystal field parameters. In other words, the 
Mj — ±4 and Mj — ±5 carry, respectively, about 90% 
of the weight of the ground state and first excited state 
doublet wavefunctions. 
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Ill) Calculation of the magnetic susceptibility of 
Tb2Ti207 and (Tbo.o2Yo.98)2Ti207. 

The d.c. magnetic susceptibility is by Van Vleck's for- 
mula [tiIJt^ using the eigenvectors determined in the pre- 
vious step. 



E 



2e, 



(2) 



where N is the number of moles of Tb'^^ , /ib is the Bohr 
magneton, fee is the Boltzmann constant, is the en- 
ergy of the i*^ level. Besides, 



.(1) - 



^(2) _ ^ mL + 2S\i;,)y 



The results as a function of the temperature are given in 
Table 3. The values for (/x(r)) are very similar for the 
Tb2Ti207 and (Tbo.o2Yo.98)2Ti207 compounds. There- 
fore any large experimentally stated difference between 
the two compounds at low temperature cannot be ac- 
counted for by a difference in the individual character- 
istics of Tb'^^ in the dense Tb2Ti207 and the dilute 

(Tbo.o2Yo.98)2Ti2 07. 



2 



TABLE 1 

Empirical CFP's deduced by the covalo-electrostatic 
model for= Tb2Ti207 and (Tbo.o2Yo.98)2Ti207, and "ex- 
perimental" CFP's for Tb2Ti207. The latter are ob 
tained using two approaches, (a) refers to the value of B5 
obtained by transposing to Tb^+ the value of Bq deter- 
mined from spectroscopic data of Eu2Ti207. (b) refers 
to the values _B^'-* obtained by transposing to Tb'^"'" the 
values of B^-^ determined from inelastic neutron data 
on Eu2Ti2 O7. All values in K. (PCEM = point charge 
electrostatic model.) 



Tb2Ti2 07 


Bl 

PCEM 471 
Cov. 610 


Bt 
708 
1599 


B\ 
-187 
-227 


Bl 
1261 


Bl 
314 




Total 1081 


2307 


-414 


1261 


314 


(Tbo.02Yo.98)2Ti2 07 


PCEM 407 
Cov. 609 


731 
1711 


-210 
-288 


1324 


389 




Total 1016 


2442 


-498 


1324 


389 


"Experimental" 


622'* 


3691'' 


-1698'' 




1665'' 



Bl 



482 
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TABLE 2 



Lowest energy levels (in K), irreducible representations in 
(Tbo.o2Yo.98)2Ti207 and Tb2Ti207, and leading compo- 
sitions of the corresponding eigenvectors. Lines labelled 
(a) and (b) are for the predicted projections for the dilute 
Tbo.o2Yo.98)2Ti207 and dense Tb2Ti207, respectively, 
using the CFP's listed in Table 1. Lines labelled (x3) 
arc for the predicted projections for the dense Tb2Ti207 
using the "experimental" CFP's listed in Table 1. 

(E) En , +/.4 ,+/-! 7f6,-/+5 ^Ge{l), +/-4 ^Ge{3) , +/-4 

(a) 0. -0.97 +/-0.06 -/+0.06 +0.15 -0.14 

(b) 0. -0.97 +/-0.G5 -/+0.G6 +0.15 -0.14 

(c) 0. -0.95 -I-/-0.13 -/4-0.13 +0.14 -0.13 



5 ^Fe ,+/-2 ^F6,-/+4 ^ 
+/-0.13 -/+0.07 
+/-0.11 -/+0.06 
+/-0.28 -/+0.14 



6(1), +/-5 5G6(3) , +/-5 
+0.14 -0.13 
+0.14 -0.14 
+0.14 -0.13 



(E) En 

(a) 18.7 

(b) 15.9 

(c) 21.6 

(A2) En 

(a) 86.5 

(b) 87.9 

(c) 85.0 

(Al) En 

(a) 121.1 

(b) 118.2 

(c) 119.6 



Fe , +/- 
-0.96 
-0.97 
-0.92 

'F6,-3 
+0.67 
+0.67 
+0.65 

'F6,-3 
+0.66 
+0.66 
+0.64 



^F6,3 
+0.67 
+0.67 
+0.65 

'F6,3 
-0.66 
-0.66 
-0.64 



^F6,-6 
+0.16 
+0.16 
+0.22 

'F6,-6 
+0.21 
+0.20 
+0.25 



'F6,6 
-0.16 
-0.16 
-0.22 

'F6,6 
+0.21 
+0.20 
+0.25 



'G6(l),-3 
-0.10 
-0.10 
-0.10 

'F6,0 
+0.07 
+0.06 
+0.11 



'G6(3),3 
-0.10 
-0.10 
-0.10 



TABLE 3 



Ionic magnetic moment (in Bohr magnetons) and inverse 
molar magnetic susceptibility as a function of tempera- 
ture. 



T (K) {/J,) (/Ub) 1/x (mole/emu) 1/x (mole/emu) 





Tbo.02Yo.98Ti207 


Tbo.02Yo.98Ti207 


Tb2Ti207 


1 


5.90 


5.8 


0.11 


5 


6.84 


21.4 


0.41 


10 


7.58 


34.8 


0.67 


15 


7.97 


47.3 


0.92 


20 


8.20 


59.5 


1.17 


30 


8.47 


83.5 


1.65 


300 


9.22 


706 


14.09 
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